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ABSTRACT 


An incompressible Jeffrey fluid heated from below in a porous medium, stability is taken into consideration, as well as how 
rotation and suspended particles may impact it. A normal mode analysis method has been used to create and quantitatively 
solve the dispersion relation. While the suspended particles are shown to destabilize stationary convection, rotation is 
found to assist stabilize the system. It is discovered that, depending on the situation, the medium permeability and the 
Jeffrey parameter can either stabilize or destabilize the system. The effects of rotation, suspended particles, Jeffrey 


parameter and medium permeability have all been depicted in graphs. 
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1 INTRODUCTION 


Convection is considered in porous medium because it has wide range of applications in fluid flow and heat transfer as 
well as heat-exchanger, oil recovery, construction materials, cooling of electronics, minimising pollutant generation and 
medicinal treatments8. Lapwood6has researched thermal convection of fluid in porous medium and Nield and Bejan 7 


have written a book on the subject. 


Rotation results in the introduction of several new elements into the problem. In rotational fluid dynamics, 
vorticity is related to a number of findings. The thermal instability of rotating fluid layers heated from below has been 
discussed by Chandrasekhar2,who also demonstrated that rotation avoids the development of instability and shows 


asymptotic behaviour. In 13,15,18,the effect of rotation on thermal instability was investigated. 


Scanlon and Segell4looked at how suspended particles affected the commencement of Bénard convection and 
discovered that coarse particles decreased the exponential growth rate of unstable disturbances while small particles 
increased it. The impact of suspended particles has been researched by 3 for more realistic boundary conditions.By heating 
fluid from below,1lexamined the impact of rotation and suspended particles in a porous media. Numerous additional 


researchers, including 1,4,12,16investigated convection with suspended particles. 


Non-Newtonian fluid has a wide range of uses in the geophysical, chemical and biological sciences. Non- 
Newtonian fluids have a wide range of industrial and technical uses, which has increased interest in their research. One of 
the simplest non-Newtonian fluid models is the Jeffrey non-Newtonian fluid, which has a time derivative rather than a 


convective derivative. Jeffrey Sused below-surface heating to study the fluid layer's stability. Jeffrey parameter has a 
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stabilising influence on stationary convection, according to recent research10on the effect of Jeffrey nanofluid in a porous 
media and 9on the same problem taking the effect of rotation into account. The start of Jeffrey fluid in a porous heat- 


generating layer has been studied in 17. 


We are interested in researching the impact of rotation and suspended particles on the stability of Jeffrey fluid 
heated from below in porous media in the current work in light of the significance of the many applications 


described above. 
2 Mathematical Model 


Consider a Jeffrey fluid confined between two parallel horizontal planes z = 0 and z = d saturated by layer of porous 


medium subjected to uniform rotation with angular velocity 2(0,0,Q) and gravity g force acting vertically. The physical 


system is heated from below such that uniform temperature 6 (= fea azl) gradient is maintained. 


g(0,0,—g) 


TrEeT TT ETT 


HEATED FROM BELOW 


ToT > 1) 


Figure 1: Physical Configuration. 


The relevant equation of state, equation of continuity, equations of motion and equation of energy for Jeffrey fluid 


with suspended particles under rotation in porous medium under Boussinesq approximation are: 


P = po — a(T — Ty) (1) 
Viqg=0 (2) 
A(2+4@.0)q) =-0p + 9(1+%)-— Ag +2 xm +2 (a -a) (3) 


where q, is the velocity of particles,qis the velocity of fluid,eis the porosity,pis the pressure,k,is the medium 
permeability, A3is the Jeffrey parameter,N is the number of densities of the particles, K = 67a is the constant given by 
Stokes drag formula, wherey is the coefficient of viscosity anda is the radius of particle.Due to the presence of suspended 
particles, there added an extra force term proportional to the velocity difference between particles and fluid in momentum 


equation (3). 


The equation of energy with suspended particles is given by 


a 
(e(pC), + (1 —£)(0C), 2 + (oC) -(q.V)T + mNG, (es —+ (@.v))T = Km V?T 


Using Boussinesq approximation 


me ot (q. v))T = «KV°T (4) 
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where,F = (« +(1-e) =. is constant, 
PoCr 


km 
Pole 


K= is the thermal diffusivity,k,, is the coefficient of thermal conductivity,(pC), is the heat capacity of 
solid,poCy is the heat capacity of fluid,mis the mass of the particles. 


The equations of motion for particles are 
a 
amN (5+ (dp-V)) dp = KN(Q — Gp) (5) 
The equation of continuity for particles is 
aN 
5, + V.(Ndp) = 0 (6) 


There must be an additional force term in the equation of motion for the particles that is same in magnitude but 
opposite in sign because the force the fluid exerts on the particles is equal to and opposite from the force the particles exert 
on the fluid. The particles' buoyancy force is disregarded. Assumedly, the distance between the particles is considerably 
greater than their combined diameters. The equations describing the movements of particles (such as equation 5) have been 


written using these presumptions. 
The initial state of the system (there is no motion) is 
q = (0,0,0), gp = (0,0,0),T — To = —62,p = po(1 + a@Bz), No = constant 
Let q(u, v,w), dp(L1r, S), 6p denotes the perturbations of velocity, particle velocityand pressure respectively. 
Let the change in temperature distribution be T' = Ty — Bz + 0 
The change in density after perturbation@ in temperature is 
6p = —@po0 


Using these perturbations in equations (2),(3) and (4) to (6) and neglecting the terms with high powers and 


products of perturbations, the resulting linearized perturbed equations are 


Vig =0 a 
10q_ __v 4 el a 
Sa ee eae (axa) Po (a — 4) ne 


where v = “is the kinematic viscosity. 


(F + eb) = B(w + sb) + xV?0 (9) 
where b = mn 
olf 
m 0 
(22+ 1) Gp =4 (10 
Vqp = 0 (11 


Now, eliminating q, from equations of motion for fluid (8) using equations of motion for particles (10), weget 
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1(m 44) 2a _ (may y)f_ 4 = mn 4 
ele ae) ge gL 5 OO ee at | oe (2) 
In cartesian form, equation (7),(9),(12) can be written as 
Ou Ov Ow 
re os a 0 (13) 
1 _ (ma 1 d6p v 2 _ ™mN Ou 
aaa at * 1)% at G att 1) | Po Ox kj fiwag) & ef E (v0)| Epo Ot Ge) 
1/ma dv (ma 1 06p v 2 __ ™N Ov 
AG att 1) at ¢ ant 1) | me tee (un)| EPo Ot (15) 
1 _ (ma 1 06p v = ™mN OW 
(28 att 1) 7 ¢ att 1) | po az kts)” - ag6| Epo Ot (16) 
Equation (9) after eliminating s becomes 
mo te] _ 2 _ mo 
(22 +1) ((F + eb) 2-Kv?)0 =B(=24+1+b)w (17) 
Operating equation (14) by a and equation (15) by < and subtracting, weget 
1/ma oF _ (ma 22 0w mN OF 
ag att 1) at G ant 1) | € Oz ky (1+A3) | Epo Ot (18) 
where € = ~ a = is the z component of vorticity. 


Further, Operating equation (14) by = and equation (15) by 5 > adding and using equation of continuity in cartesian form 


(i.e., (13)), weget 
ee a Get Meee Ot (19) 


Now, operating equation (19) by 2 and using equation (16) to eliminate dp , weget 


ater = 5. V*w) = (Ga a 1)| "aE mw : - ae (= = =) 6| a Z a oe) 


Boundary conditions at free surfaces: At free surface tangential stresses are zero. 
w=0,0=0atz=0andz=d; (20) 


azw 
Oz? 


= 0,¢ = 0 ona free surface 
3.Normal Mode Analysis 
We attributea dependency on x, y and t of the form 


exp(i(k,x +kyy) + nt) , where, k, and k, are the wave numbers along the x — axisand y — axis respectively, k? = 


ke + he is the resultant wave number and 7n is the complex constant known as growth rate. 
We consider the perturbations ofw, 0,¢ having the form’ 


[w, 6,7] = [W(z), 0(z), E(z)lexp(i(k,.x + kyy) + nt) (21) 
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Using equation (22) in equations (17), (18) and (20), weget 


2(BR+1) (S-w)w = (BL41) pag (S- ew 2 gato 


ee ae ie eee aa 


(FS+ 1)(@ + 26)n~«(S- ))e =p(ES+1+b)W 


It is convenient to discuss equation (23)-(25) in non-dimensional variables. 


nd? m Tv mN 


—,T,=—,M= 


=x'd,y =y'd,z=2'd,a=kd,o=—,t= 
x=xdy=yd,z=z'da ae ae 


a2’ Pe ee eal 
0 
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(22) 


(23) 


(24) 


where the non-dimensional parameters are: PR, = zis the medium permeability, p, = ~ is the Prandtl number. 


[E(1 + M ) + | (D2 a®)w + 222 pg 4 eee d? 


TyO+1 


[E(a +—“)+ ater - a 


2 
(D2 — a? — Eyp,0)0 = oe Ww 


K Tyo+1 


0=0 


Eliminating © and € from equation (26) using equation (27) and (28), weget 


He + 4) + mon (D? = a? = Eyp,a)(D? - a?)W + 5 (D? — a? - Eyp,o)D?W = Ra a(act’) He + 


To+1 


Where T = a is the Taylor number and R = ace d* is the Rayleigh number. 


Boundary conditions (21)using non-dimensional variables becomes 


W =0,DE =0atz=O0andz=1 


and D*W = 0,0 = 0 ona free surface 


The proper solution of W at lowest characteristic satisfying the boundary conditions (30) must be 


W =W,sinzz 
Where W, is constant. 


Using this solution in equation (29), weget 


Ra? (222) = [21 +) + | (ne? +a? + Epo) (0? + a?) 


Tyo+1 é Ty0+1 
5 nT (1? + a? + Epo) 
eG Gea) ea) 
Further eliminating zz by letting 


a? o R T 
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Ra= Ig aga Py Tz = UT, 


(25) 


(26) 


(27) 


(28) 


(29) 


(31) 


(32) 
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Weget 


Ryx = (22) (26 +) tap) ht 2 + pro) (1 + x) + ee | (30) 


it20,+B € T204+1 Py(1+A3) ( - ( oereTe PiGtaa) 


4. Stationary Convection 


Put o = 0, for stationary convection equation (33) reduces to 


Ra 


_ (1+x) 14x moh) 
~ xB LPy(1+43) e2 


(31) 


Which expresses the modified Rayleigh number as a function of the dimensionless wave number, Taylor numberT,, 


medium permeabilityP,, porositye, Jeffrey parameterA3, suspended particle parameterB. 


ORq ORq ORq ORq 
OB’ AP,’ ATq’ OAz 


To study the effect of B, P,,T,,A3and €, we examine the behaviour of and “ + analytically. 


Equation (34) yields 


ORqg _—(Atx)[ 1tx TaP(1+43) 

OB xB? LP(1+A3) E2 | (32) 
Which is negative so the suspended particles have destabilizing effect on the physical system. 
From equation (34), we get 

ORq _ (+x) = 1+x Tg(1+A3) 

aP, xB P12(1+A3) e2 | (33) 

: oe Sas ts 1+x Ta(1+A3) ae 1+x 
Thus, the medium permeability has stabilizing effect when GHD eee and destabilizing effect when reer 
Tq(1+43) 
ez 

From equation (34), we also get' 

ORq _ (14x) [Py (1+A3) 

aT. «XB e2 | (34) 
This implies that rotation has a stabilizing effect on the system. 
It is evident from equation (34) that 

ORq _ (+x) 1+x TaP1 

aa3..—~—t«X#@’ Py (1+A3)2 + e2 | (35) 

atts we tx TaP1 ape we 14K TaP1 

Thus, the Jeffrey parameter has stabilizing effect on system if Patan? ar and destabilizing effect if pus? 
Equation (34) yields 

ORq _ _2(1+x) [TaPy(1+As) 

Qe xB 3 | (36) 


Porosity has destabilizing effect on the physical system. 
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5. Principle of Exchange of Stability 


Multiplying equation (26) by W* and integrating over the range of z (i.e., from z = 0 to z = 1) and making use of 


equation (27) and (28), weget 


{f+ )+— Jn -¢ [40+ +2 gaa (2241), + o*p,Eyla] = 0 (37) 


E Tyo +1 Pm(1+A3). T10*+1 Pm(1+A3) vB T10*+B 


where, 
1, = f, (Dwi? + a*|W|?)dz 
1izy2 
I, = fy \§|?dz 
Is = f) (\DO|? + a?|@|?)dz 
le@y2 
I, = f, |Ol’dz 
The integrals are all positive definite. 


Puttingo = io;in equation (40), and equating imaginary parts weget 


1 M 2 gaa*x [ t1(B-1) (ty0))°+B = 
oi tt (1 = TOEND li, + d*l] + vB feaee 3" (10;)2+B2 opr Exl4]} = (38) 
0; may not always be zero. i.e., 
1 M 2 gaa*x [ t,(B-1) (t10;)?+B = 
é (1 - aaa) els vB Freer: 3°" (t10))2+B2 oiprEyLa| =e 


Which gives the possibility that 0; # 0 


As a result, given the existence of rotation, suspended particles and the Jeffrey parameter, the exchange of 


stabilities concept may not be applicable to this issue. 
6. Numerical Results 


The variation of thermal Rayleigh number for stationary case with suspended particles, medium permeability, Taylor 


number, Jeffrey parameter and porosity for fixed wave numbers have plotted using equation (34). 


Figure 2 shows the variation of R,for stationary convection with suspended particles B for different values of 
wave numberx = 0.1,0.5,0.9 and fixed values ofP, = 0.2,T, = 100,A3 = 0.6 and ¢ = 0.5. As the values of Bincreases 


the graph shows downward slope, thereby destabilizes the stationary convection. 
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Figure 2: Variation of R, with B for three values of wave numberx = 0.1,0.5,0.9 
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Figure 3: Variation of R, with P, for three values of x = 0.1,0.5,0.9. 


In figure 3, R, is plotted against the medium permeability P, for different values of wave number x = 0.1, 0.5, 0.9 


and fixed values of B = 10,T, = 100,A3 = 0.6,¢ = 0.5 and has both destabilizing and stabilizing effect. 
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Figure 4: Variation of R, with T, for three values of x = 0.1,0.5,0.9. 
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In figure 4, R, is plotted against the Taylor number T, for different values of wave number x = 0.1, 0.5, 0.9 and 


fixed values of B = 10, P, = 0.2,A3 = 0.6,¢ = 0.5. It is clear from the graph that effect of rotation has stabilizing effect 


on stationary convection. 


Figure 5: Variation of R, with A; for three values of x = 0.1,0.5,0.9. 


In figure 5, the variation of R, for stationary convection with Jeffrey parameter A3 for different values of wave 


numberx = 0.1,0.5,0.9. and fixed values ofB =10,P,=0.2,T, =100 and e¢=0.5. The graph shows 


destabilizing/stabilizing effect on the physical system. 


Rs 
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Figure 6: Variation of R, with € for three values of x = 0.1,0.5,0.9. 


In figure 6, Rg is plotted against the porosity,e = 0.5. It is clear from the graph that effect of porosity has 
destabilizing effect on the physical system. 


7. CONCLUSION 


In the current study, we investigated the linear stability theory on the start of the Jeffrey fluid by heating it from below in a 
porous media while taking into account the effects of rotation and suspended particles. For stationary convection, 


analytical and graphical solutions have been obtained and we draw the conclusion that 


1. Taylor number on stationary convection show an asymptotic behaviour and has a stabilizing effect on the system. 


1+x Ta(1+d3) 


2. The medium permeability has stabilizing effect when pete) 2 


and destabilizing effect when 
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14x Tg(1+A3) 
P2(1+A3) e2 


Suspended particles Band medium porosityehas destabilizing effect on the stationary convection. 


The Jeffrey parameter has stabilizing effect on stationary convection if ae < “an 
1 3 


and destabilizing effect if 


1+x TaPy 
Py (1+A3)2 B2° 


The principle of exchange of stabilities may not be valid in this problem due to the presence of rotation, 


suspended particles and Jeffrey parameter. 


Acknowledgements 


Authors would like to thank the reviewers for their valuable suggestions and comments for the improvement of quality of 


the paper. 


REFERENCES 


1. 


10. 


11. 


Aggarwal, A.K.: Effect of rotation on thermosolutal convection in a Rivlin-Ericksen fluid permeated with 


suspended particles in porous medium. Advances in Theoretical and Applied Mechanics, 3, 177-188 (2010) 
Chandrasekhar S.: Hydrodynamic and Hydromagnetic stability. Dover Publication, New York (1981) 


Chand, R., Rana, G.C., Hussein, K.: Effect of suspended particles on the onset of thermal convection in a 
nanofluid layer for more realistic boundary condition. International Journal of Fluid Mechanics Research, 42, 


375-390 (2015) https://doi.org/10.1615/InterJFluidMechRes. v42.i5.10 


Chand, R., Rana, G.C., Kumar, A., Sharma. V.: Thermal instability in a layer of nanofluid subjected to rotation 
and suspended particles. Research Journal of Science and technology, 5, 32-40 (2013) 


Jeffrey, H.: The stability of a layer of fluid heated below. Journal of Science, 2, 833-844 (1926) 


Lapwood, E.R., Convection of a fluid in a porous medium. Proceeding of the Cambridge Philosophical Society, 
44, 508-521, (1948). https://doi.org/10.1017/S030500410002452X 


Nield. D.A., Bejan. A.: Convection in porous medium. Springer, New York (2013) 


Nield, D.A., Simmons, C.T.: A brief introduction to convection in porous media. Transport in Porous Media, 130, 


237-250 (2019). https://doi.org/10.1007/s11242-018-1163-6 


Rana. G.C.: Effect of rotation on Jeffrey nanofluid flow saturated by porous medium. Journal of Applied 
Mathematics and Computational Mechanics, 20, 17-29 (202 I)https://doi.org/10.17512/,amcm.2021.3.02 


Rana, G.C., Gautam, P.K.: On the onset of thermal instability of a porous medium layer saturating a Jeffrey 
nanofluid. Engineering Transactions, 70, 123-139 (2022). https://doi.org/10.24423/EngTrans.1387.20220609 


Rana, G.C., Kumar, S.: Effect of rotation and suspended particles on the stability of an incompressible Walters’ 
(Model B) fluid heated from below under a variable gravity field in a porous medium. Engineering Transactions, 


60, 55-68 (2012) 


Impact Factor (JCC): 6.6810 NAAS Rating 3.45 


Effect of Rotation and Suspended Particles on the Stability of Jeffrey Fluid in a Porous Medium 67 


12. 


13. 


14. 


15. 


16. 


17. 


18. 


19. 


Rana, G.C., Thakur, R.C.: Combined effect of suspended particles and rotation on double-diffusive convection in 
a viscoelastic fluid saturated by a Darcy-Brinkman porous medium. The Journal of Computational Multiphase 


Flows, 5, 101-113 (2013). https://dot.org/10.1260/1757-482X.5.2.101 


Rana, G.C., Thakur, R.C.: The onset of double-diffusive convection in a layer of nanofluid under rotation. Revista 
da Engenharia Térmica, 15, 88-95 (2016). http://dx.doi.org/10.5380/reterm.v15i1.62153 


Scanlon, J.W., and Segel, L.A.: Effect of suspended particles on onset of Bénard convection. The Physics of 
Fluids, 16, 1573-1578 (1973). https://doi.org/10.1063/1.1694182 


Sharma, R.C.: Effect of rotation on thermal instability of a viscoelastic fluid. Acta Physica Academiae 
Scientiarum Hungaricae, 40, 11-17 (1976). https://doi.org/10.1007/BF03157148 


Sharma, V., Gupta, S.: Thermal convection of micropolar fluid in the presence of suspended particles in rotation. 


Archives of Mechanics, 60, 403-419 (2008)https://doi.org/10.5958/2349-2988.2017.00037.7 


Yadav, D.: Effect of electric field on the onset of Jeffrey fluid convection in heat-generating porous medium layer. 
Pramana-J Phys, 96, 19 (2022). https://doi.org/10.1007/s 12043-021-02242-6 


Yadav, D., Bhargava, R.: Thermal instability of rotating nanofluid layer. International Journal of Engineering 
Science, 49, 1171-1184 (2011).https://doi.org/10.1016/j.ijengsci.2011.07.002 


Wooding, R.A., Rayleigh instability of a thermal boundary layer in flow through a porous medium. Journal of 
Fluid Mechanics, 9, 183-192 (1960). https://doi.org/10.1017/S0022112060001031 


www. laset.us editor@iaset.us 


